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Abstract. In the Hopf algebra approach of Connes and Kreimer on renormalization of 
quantum field theory, the renormalization process is views as a special case of the Algebraic 
Birkhoff Decomposition. We give a differential algebra variation of this decomposition and 
apply this to the study of multiple zeta values. 

1. Introduction 

This paper applies the renormalization method in quantum field theory to the study of 
multiple zeta values when the defining sums of the multiple zeta values are divergent, with 
an emphasis on the differential structure underlying the renormalization process. 

As we will see later, such divergence of the infinite sums cannot be cured by traditional 
mathematical methods, such as analytic continuation as in the case of one variable (Rie- 
mann) zeta function. On the other hand, theoretical physicists have dealt with similar 
divergencies in quantum field theory for several decades with such a great success that the 
related renormalization process is regarded as one of the greatest achievements in mod- 
ern physics. Nevertheless mathematicians have been skeptical about the soundness of the 
mathematical foundation of the renormalization process. This situation is changed by the 
recent seminal work of Connes and Kreimer [T2J [13] which puts the renormalization process 
in a more mathematical framework of Algebraic Birkhoff Decomposition. Such a framework 
not only reveals the mathematical structure underlying the physics process, it also makes it 
possible for this physics process to be adapted to treat other apparently unrelated problems 
in mathematics, in particular the divergence of multiple zeta values. 

This is the goal of this and the companying paper [22] which complete each other. This 
paper consists of three sections. Section [2] motivates our renormalization approach through 
examples and special cases. In Section [31 we obtain the differential variation of the Algebraic 
Birkhoff Decomposition of Connes and Kreimer [12], [13] . We also extend the well-known 
quasi-shuffle Hopf algebra of Hoffman [32] in two directions. In one direction we consider 
such an algebra that is generated by a module instead of by a locally finite set. In the other 
direction, we consider quasi-shuffle Hopf algebras in the context of differential Hopf algebras 
instead of Hopf algebras. This context allows us to apply the Differential Algebraic Birkhoff 
Decomposition in Section H] to study differential properties of renormalized multiple zeta 
values. 

This paper highlights the differential aspect of the renormalization process. As is well- 
known, for the convergent multiple zeta values, the interaction between the quasi-shuffle 
product from their summation representations and the shuffle product from their integration 
representations plays a major role in their study. In the renormalized approach of divergent 
multiple zeta values, the integration representation is not available. On the other hand, 
differentiation is essential in the study of renormalized multiple zeta values. We interpret 
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this in the Algebraic Birkhoff Decomposition and give an application on generating functions 
of renormalized multiple zeta values. 

Parts of this paper, such as the examples in Section [2] and the general construction of 
quasi-shuffle Hopf algebras in Section [31 were included in an early version of [29] , but was 
taken out of the final version of that paper to limit the size of the paper. Some readers 
have found the examples helpful in understanding the renormalization process in general. 
Further the general construction of quasi-shuffle algebra in Section [3] has been referred to 
or related to in a number of papers, such as that of D. Manchon and S. Payche [37] on 
renormalization of multiple zeta values, of J. Zhao [37] on renormalization of multiple q- 
zeta values and of F. Menous [30] on Birkhoff decomposition. In particular, for the last two 
references, the case treated in [29] is not adequate and the general construction is needed. 
Also considering its application in later part of this paper in the differential context, we 
have decided to make the details of the construction available. 

Acknowledgements. Both authors thank the Max-Planck Institute for Mathematics in 
Bonn for the stimulating environment. The first named author thanks the NSF grant DMS 
0505445 for support. Thanks also go to Robert Sczech, especially for suggesting the term 
algebraic continuation, and to Herbert Gangl, especially for suggesting the bar notation. 



where si, • • ■ , s& are positive integers with s\ > 1. Since the papers of Hoffman [30] and 
Zagier [3S] in the early 1990s, their study have attracted interests from several areas of 
mathematics and physics [7J El EU [221 [231 E21 E3] , including number theory, combinatorics, 
algebraic geometry and mathematical physics. 

In order to study of the multiple variable function C( s i; ■ • • > s k) a ^ integers si, • • • , s& 
where the defining sum ([1]) is divergent, one first tries to use the analytic continuation, 
as in the case of the one variable case of the Riemann zeta function. Such an analytic 
continuation is achieved recently in [H [31 [381 EHl EE]. Unfortunately, unlike in the one 
variable case, the multiple zeta function in Eq. (pP) is still undefined at most non-positive 
integers even with the analytic continuation. Nevertheless, possible definitions of multiple 
zeta functions at certain non-positive integers were proposed in [H |3] by making use of the 
analytic continuation. Let us briefly recall these previous progresses before introducing our 
approach by renormalization. 



2. Renormalization of MZVs: motivation and examples 



Multiple zeta values (MZVs) are defined to be the convergent sums 



(1) 




2.1. Earlier approach by analytic continuation. Analytic continuation of £(si,S2) 
has been considered as early as 1949 by Atkinson [5] with applications to the study of the 
asymptotic behavior of the "mean values" of zeta-f unction. 
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Through the more recent work of Zhao [46] and Akiyama-Egami-Tanigawa [I], we know 
that C( s i; ' ' ' > s k) can be meromorphically continued to C fc with singularities on the subva- 
rieties 

s 1 = 1; 
si + S2 = 2, 1, 0, — 2, — 4, • • • ; and 

3 

e Z< 3 - (3<j<k). 

i=l 

We also have some control of the behavior near the singularities. For example, near (0, 0), 

j-t v 5si + 4s 2 

CK, S 2 ) = —, r + R2K8X, S 2 ) 

12(si + s 2 J 

where -R 2 (si, s 2 ) is an entire function near (0, 0) with -R 2 (0, 0) = 0. 

In [U [3], several definitions were proposed for multiple zeta functions at (si, - • • , s^) 
where Sj are all non-positive. Some of the definitions of C( s i ; ' ' ' , s k) ar e 

(2) lim ■ • • lim ((r 1 , ■ ■ ■ ,r k ), lim • • • lim C(n, • ■ ■ , r k ), lira C(si + r, ■ ■ ■ , s k + r). 

Naturally they give different values. In the case of C( s i> ^2) at (si, s 2 ) = (0, 0), the proposed 
values are 5/12, 1/3,3/8 respectively according to the above three definitions. In fact, by 
letting (7*1, r 2 ) approach to (0, 0) along different paths, one can get any value, as well as the 
infinity, to be the limit. Even though some good properties of the variously defined non- 
positive MZVs were obtained in these studies, they fell short of the analogous properties of 
the positive MZVs, especially the double shuffle relations. 

2.2. An illustration of the renormalization approach. Our approach to define MZVs 
where the sum (TjQ) is divergent is adapted from a renormalization procedure (dimensional 
regularization plus minimal subtraction) in quantum field theory (QFT). This definition 
coincides with the usual definition of MZVs when they are defined. Our extended MZVs 
also satisfy the quasi-shufffe relation. So our approach is an algebraic continuation in 
the sense that we extend the definition of MZVs that preserves the quasi-shuffle relation. 

The QFT renormalization procedure was recently put into a more mathematical frame- 
work through the work of Connes and Kreimer [121 C3] an d was thus made it possible to 
be applied to problems in mathematics. For our purpose, the dimensional regularization 
of Feynman integrals is replaced by a regularization (or deformation) of infinite series that 
has occurred in the study of Todd classes for toric varieties [H]. 

Here we illustrate this method by some special cases with the general cases given in §H] 
and in [25] . First recall the following generating series of Bernoulli numbers that goes back 
to Euler. 

(3) ?V5> J 

fc>0 

It can be easily rewritten as 

(4) — = ---^ = -- + VC(-A/ 



e 
AT 



k>0 

since B = 1 and ((-k) = (-l) k ^± for k > 0. 



k\ 



4 LI GUO AND BIN ZHANG 

Now consider 



n>l 

regarded as a deformation or "regularization" of the series defining the Riemann zeta func- 
tion £(s) = J2 n >i The regularized series converges for any integer s when Re(e) < 0. 
In particular, 

1 — e e 

and Eq. (jlj) gives the Laurent series expansion of the regularized sum Z(0; e) = Y^n>i eUe a ^ 
e = 0. For a Laurent series f(e) = ^2 n>N a n e n where N G Z, we denote the power series 
part (or the finite part) Xm>o a « £ri of the Laurent series by fp(/), a notation borrowed 
from [36] which, according to the authors, can be traced back to Hadamard. We tentatively 
call fp(/)| e =o the renormalized value of f(e) (see Section HJ] for the general case). Then we 
have 

f p(E e " £ ) =c(°)- 

n>l 

So the renormalized value of Z(0;e) = J2n>i en£ * s C(0)- Similarly, to evaluate £(— A;) for 
an integer k > 1, consider the regularized sum 

z(-*., e ) = 5>v»=* (T ^) 

n>l 

which converges uniformly on any compact subset in Re(e) < 0. So its Laurent series 
expansion at e = is obtained by termwise differentiating Eq. (jlj), yielding 

OO j 

(5) Z(-k; e) = (-l)- k -Hky. e'"' 1 + £ C(~* - 

3=0 J ' 

We then have 

f P (5> fe e-)| £=o = C(-A;). 

n>l 

Thus the renormalization method does give the correct Riemann zeta values at non- 
positive integers. We next extend this to multiple zeta functions and "evaluate" £(0, 0), for 
example, by consider the regularized sum 

Z(0Ae)= £ e"^ = i 5 e£i f e2£ . 

ni>'«2>0 

Naively taking the finite part as in the one variable case, we find 
fp( £ e — e — )|_ o = fp (_l_ _ | C(0 )I + ( - | C (-1) + C(0) 2 ) + o(s))|_ = 11/24. 

ni>'«2>0 

This does not agree with any of the previously proposed values of £(0, 0) in Eq. (j2J). Further, 
this value does not satisfy the well-known quasi-shufne (stuffle) relation: 

C(o)C(o)^2C(o,o) + c(o) 

since the left hand side is 1/4 and the right hand side is 5/12. Recall the well-known 
quasi-shuffle relation 

((s)((s) = 2(( S ,s) + ((2s) 
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for any integer s > 2. 

To improve this situation and obtain a more suitable definition of £(0, 0), we recall that 
a key principle in a renormalization procedure of QFT is that if a divergent Feynman 
integral contains a component integral that is already divergent, then the divergency of the 
component integral should be removed before removing the divergency of the integral itself. 
For our example of £(0, 0), the regularized sum 

e ni£ e n2e = Y e?l2£ Y e " X * 

ni>TL2>0 ri2>l m>n2+l 

has a component sum X^n >n 2 +i eUl£ ^ na ^ * s a l rea dy divergent when e goes to 0. By the 
renormalized process adopt to our case (see Eq. fl32l) and Remark 14.71) we found that the 
renormalized value should be defined by 



ni>ri2>0 n2>0 ni>0 ni>0 



e=0 



sub divergence 

fp(^ - §C(0)| + (-|C(-1) + C(0) 2 + o(e)) - (I - C M _ C (_i) + o{£ 



£=0 



= -|c(-i) + C(o) 2 = |- 

This value indeed satisfies the quasi-shuffle relation C(0)C(0) = 2C(0, 0) + C(0). 

This renormalization process will be systematically carried out for all C( s i ; ' • • ; s k) in §111 
after the general setup in §0 Remark 14.71 shows how the above example follows from the 
general case. 

3. Differential Birkhoff decomposition and quasi-shuffle Hopf algebras 
We formulate a general setup for our later applications to MZVs. 

3.1. Differential Birkhoff decomposition. We review and extend the algebraic frame- 
work of Connes and Kreimer for renormalization of perturbative quantum field theory. For 
further details of physics applications, see [T2| [T3| f!4"l fl8| fl9l l21j 135] . 

Let k be a unitary commutative ring that we usually take to be 1 or a subring of C. In the 
following an algebra means a unitary k-algebra unless otherwise specified. A connected 
filtered Hopf algebra is a Hopf algebra (H, A) with k-submodules H n , n > 0, of H such 
that 

(1) H n C H n+1 ; 

(2) U n > H n = H; 

(3) H p H q C H p+q ; 

(4) A(H n )C® p+q=n H p ®H q . 

(5) (connectedness) H = k. 

A Rota— Baxter algebra 0, H21 H3] of weight A is a pair (R, P) where R is an algebra and 
P : R — > R is a linear operator such that 

(6) P(x)P(y) = P(xP(y)) + P(P(x)y) + XP(xy), 

for any x, y G R. Often 6* = — A is used, especially in the physics literature. 
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Now let if be a connected filtered Hopf algebra and let (R, P) be a commutative Rota- 
Baxter algebra. Consider the algebra 31 := Hom(if, R) of linear maps from if to R where 
the product is given by the convolution 

f-kg(x) := 53/(x(i))y(z( 3 )). 

(x) 

Here we have used the Sweedler's notation A(x) = J2( x ) x (i)® x {2)- Then the operator 

? : Hom(ff, R) -> Hom(ff, 4/hPo/, 
is a Rota-Baxter operator on 3?. 

Definition 3.1. (1) A differential algebra is a pair (R,d) where R is an algebra 
and d is a differential operator, that is, such that d(xy) = d(x)y + xd(y) for all 
x,y G R. A differential algebra homomorphism / : (Ri,d\) — > (^2,^2) between two 
differential algebras (i?i, di) and (i?2, ^2) is an algebra homomorphism / : R x — > R 2 
such that / o di = d 2 o f. 

(2) A differential Hopf algebra is a pair (H, d) where if is a Hopf algebra and 
d : H — > H is a differential operator such that 

(7) A(d(x)) = ^ i d ( x (i))® x (2) + x (i)®d(x {2 ))) ■ 

Or) 

(3) A differential Rota-Baxter algebra is a triple (R, P, d) where (R, P) is a Rota- 
Baxter algebra and d : R — > ff is a differential operator such that P o d = d o P. 

Theorem 3.2. Lei H be a commutative connected filtered Hopf algebra. Let (R,P) be a 
Rota-Baxter algebra of weight —1. Let <fi : H — > R be an algebra homomorphism. 

(1) (Algebraic Birkhoff Decomposition) There are algebra homomorphisms <p_ : 
H -> k + P(i?) and 0+ : if -> k + (id R - P){R) such that 

<p = <j> + . 

Here (f)*}' 1 ^ is the inverse of </>_ with respect to the convolution product. Further, 

(8) 0_(x) = -P(<f>{x) + J2<f>4 x ')<f>( x ")) 

(x) 

and 

(9) <f> + (x) = (id - P) + M^OO) ■ 

(x) 

Here we have used the notation A(x) = x®l + l®x + x> ® x " ■ 

(2) If P 2 = P, then the decomposition in is unique. 

(3) (Differential Algebraic Birkhoff Decomposition) If in addition (if, d) is a 
differential Hopf algebra, (ff, P, d) is a commutative differential Rota-Baxter algebra, 
and is a differential algebra homomorphism, then 0_ and <p + are also differential 
algebra homomorphisms. 
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Proof. For the proof of <^), see [35l Theorem II. 5.1]. For the proof of ([2]), see [201 Theo- 
rem 3.7]. So we just need to prove ([3]). We only verify that </>_ is a differential algebra 
homomorphism, the proof for <j) + being the same. For this we use induction on n to prove 

(10) do<p„( x ) = </>_ od(x) 
for x G H n , n > 0. 

For x G Hq, we have x G k. Since 0_ is an algebra homomorphism by ((H), we have 
<Mz) = So ^-(z)) = = (/>-{d{x)). 

Assume Eq. (TTOT) has been verified for n < k and consider x G Hk+i- Applying d to 
Eq. ([8]), we have 

0(M*O) = 9(-p(0(x) + ^0„(x')0(x"))) 

(x) 

= -P(9(0(x)) + J](<90_(«0Or") + 0_(*W(OO))) 

(x) 

= -Pfa(d(s)) + ^(0„(rf(x'))0(x") + <M*')0(d(*"))). 

(*) 

Here we have used the induction hypothesis in the last step. On the other hand, 

4>-(d(x)) = -P(0(d(*)) + 

(d(x)) 

Further 

A(tZ(s)) = d(x)®l + l®d(ar) + d(x)'®d(x)" 

(«*(*)) 

and 

d(A(x)) = + l®z + ^ = + l®d(x) + ^ {d(x')®x" + x'®d{x")) . 

(x) (x) 

Since i7 is a differential Hopf algebra with the operator d, by Eq. (CO), we have 

^ d{x)'®d{x)" = (d(x')®x" + x'®d(x")) 
(«*(*)) (*) 
and then applying 0_ * to it, we get 

(d(x)) (x) 

This completes the induction for Eq. (I10p . □ 

3.2. Mixable shuffle algebras. Let A be a unitary or nonunitary k-algebra. Consider 
the tensor product algebra 

(11) <K A :=T(A) = @ n > A* n 

with the convention that A^° = k. Here we have used | instead of <g> to denote the tensor 
product of A with itself since ® will be used for the coproduct introduced below. Let 1 be 
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the product in T(A), so for a E A* m and b G A 1 



we have 



(12) 



a |b 



a|b G A* m+n , 
ab G A* n , 
ab G Al m , 

ab G k, 



if m > 0, n > 0, 
if m = 0, n > 0, 
if m > 0, n = 0, 
if m — n — 0. 



Here the products in the second and third case are scalar product and in the fourth case is 
the product in k. Thus, 1 identifies k %A and A |k with A by the structure maps k | A — ► A 
and A |k — > A of the k- module A. 

We next equip CKa = T(A) with another product * so that, together with the deconcate- 
nation coproduct, "Ha is a connected filtered Hopf algebra. In a special case, this recovers 
the quasi-shuffle Hopf algebra of Hoffman [22] • The product * is defined by the quasi- 
shuffle product in its recursive form, and by the mixable shuffle product in its explicit 
form by one of the authors and Keigher [27]. Given two pure tensors o G A ' m and b G A ' n , 
to define a*b explicitly, recall that a shuffle of a and b is a tensor list of aj and bj without 
changing the order of the a^s and bjS. A mixable shuffle is a shuffle in which some (or none) 
of the pairs Oj §bj are merged into a>ibj. See [27] for the precise definition. Then a * b is 
the sum of mixable shuffles of a and b. It is shown in [15] that the quasi-shuffle product 
and the mixable shuffle agree with each other. We note that the mixable shuffle arises 
from the construction of free commutative Rota-Baxter algebras and has been denoted by 
o + in the literature of Rota-Baxter algebras (previously known as Baxter algebras), such 
asJlfiailSlinirjllEaESlEZlES]. We adapt the notation * to be consistent with the 
current convention on multiple zeta values. 

We give more details of the recursive definition of the product *. We only need to define 
a * b for pure tensors a G A ' m and b G A ' n , m, n > 0, and then to extend it to a product 
on "Ha by biadditivity. If m = 0, then o G k. Then we define a* b to be the scalar product. 
Similarly if n = 0. In particular, the identity 1 of k is also the identity for *. So we only 
need to define a* b when m > and n > 0. Then we have a = ai|o' and b = b\\b' for 
ai,&i G A and a' G A^ m ~ l \ b' G A^ n ~ l \ With this notation, we define a* b by induction 
onm + n. When m + n = 0, this has been done. Assuming that a* b has been define for 
< m + n < k and consider a G A^ m , b G A^ n with m + n — k + 1. If either m or n is 
zero, then again we are done. If m > and n > 0, then define 



By the induction hypothesis, the three terms on the right hand side are defined. 

Without resorting to the notation 1 , the product * is defined by the following recursion. 
First define the multiplication by A^° = k to be the scalar product. In particular, 1 is the 
identity. For any m, n > 1 and o := a\ | • • • |a m G A^ m , b := b\ | • • • \b n G A^ n , define 
a * b by induction on the sum m + n. Then m + n > 2. When m + n = 2, we have a = ai 
and b — b\. Define 



(13) 



a * b = a x | (a' * b) + &i |(a * b') + (aA) I (a' * b'). 



(14) a * b = a\ |&i + 6i |ai + aifoi. 
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Assume that a*b has been defined for m+n > k > 2 and consider a and b with m+n = k+1. 
Then m + n > 3 and so at least one of m and n is greater than 1 . Then we define 

a*b = ai|&i|--- |6 n + 6i|(a 1 *(6 2 |--- |6«)) 

(15) +(ai&i) | &2 ■ • • • l&n, when m = 1, n > 2, 
a*b = ai |((a 2 | • • • |a m ) * foi) + 6i |ai | • • • |a m 

(16) 1^2 I • • • l«m, when m > 2, n = 1, 

a * b = ai | ((a 2 1 |a m ) * (61 1 ■ • ■ |&„)) + 61 1 ((a x | |a m ) * (& 2 | • ■ ■ |& n )) 

(17) +(a 1 ,6 1 )((a 2 I |a m ) * (62 I ■ • ■ |& n )), when m,n>2. 

Here the products by * on the right hand side of each equation are well-defined by the 
induction hypothesis. 

To relate to the quasi-shuffle algebra of Hoffman and to the algebra we will use later, we 
introduce the following concept. 

Definition 3.3. A semigroup X is called a filtered semigroup if there is an increasing 
sequence X n C X, n > 1, of subsets of X such that X m X n C X m+n , m,n > 1. 

Let Ax = k X be the semi-group ring of X with coefficients in k. A special case is when 
X is locally finite, i.e., X has a grading given by the disjoint union X = ]J fc>1 X^ with 
each finite and such that X^X^ C x( m+n \ m,n > 1. 

The following theorem is a simple generalization of [521 Theorem 3.1, 3.2] and [TBI The- 
orem 2.5] where X is a locally finite set. The same proof works in our generality with the 
notation |. 

Theorem 3.4. Let A be a commutative nonunitary k-algebra. 

(1) The product * equips J~C A with the structure of a commutative unitary filtered k- 
algebra. 

(2) The product * coincides with the mixable shuffle product in [27] . 

(3) Equip *Ka with the deconcatenation coproduct 

A : [K4 — > ^^(8)3^,4) 

m— 1 

A(ai I • ■ ■ |a m ) = l®(ai | |a m ) + ^(fli I lai)®{a i+1 1 |a m ) 

i=l 

(18) +(ai 1 ■ • ■ |a m )®l 
and t/ie projection counit 

e : 'Ha — ► k 

onto t/ie direct summand k = A '° C Dt" A . T/ien IK^ zs a commutative cocommutative 
connected filtered Hopf algebra. 

Proof. (TTJ The same proof for [321 Theorem 2.1] applies. Just replace the length of a 
word by the tensor degree of a tensor. Alternatively, one can use item (J2]) which is proved 
independent of item ([T]) and use the fact that the mixable shuffle product is commutative 
and associative [27J Theorem 3.5]. 

(J2J) The proof is the same as for [T5l Theorem 2.5] since the proof there only requires that 
A is a commutative k-algebra. 
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([3]) The proof that 'Ha is a bialgebra is the same as for [321 Theorem 3.1] by considering 
tensor products in place of word concatenations. Since the result is essential for our later 
applications, we provide some details. 

The coassociativity is clear. So to prove that %a is a bialgebra, we only need to show 
that e and A are algebra homomorphisms. For e, this is clear. For A, we prove 

(19) A(o) * A(b) = A(a* b) 

for pure tensors a E A^ m and b G A ' n by induction on m + n > 0. If m + n < 1, then at 
least one of a and b is a scalar in k and Eq. ffl9|) is clear. Suppose the equation has been 
proved for < m + n < k and consider pure tensors a and b with m + n — k + 1. Then 
Eq. f fl9|) is again clear if either one of m or n is zero. So we can assume m > and n > 0. 
Then we have a = ai|a' and b = &i|b' with aiA G A and a' G A^" 1 " 1 ), b' G A ,(n ~ i:) . Let 

A(o') = a (D ® a (2)> A ( b ') = Yl b (D ® b (2) 

(a') (b') 

by Sweedler's notation. Then by Eq. ([TBI) , we have 

(20) A(o) = 1 ® o + 5^(aiio' (1) ) ® 0J2), A(b) = 1 ® b + J^(6i |b' (1) ) ® b' (2) . 

(<*') (V) 

That is 

(21) A(a) = l®a + a 1 |A(a / ), A(b) = 1 ® b + 6 1 |A(b / ). 

Here ai|A(o') means a\ is multiplied with the first tensor factors of A(a'). Thus 

A(a)*A(b) = ^((aila / (1) )*(6 1 ib , (1) ))®^( 

a (2) * b (2)) 

(a') (b') 

+ Xl( ai ' a (i)) ® ( a (2) * b ) + 5Z( & il b (i)) ® ( a * b (2)) + 1® (a* b). 

(a') (b') 

Applying Eq. ({TBI to the first and fourth terms gives 

( a iK a (i) * ( fo il b (i))) + & ii(( a il a (i)) * b (i)) + (ai&i)lKi) * b (i))) a (2) * b (2)) 

(a') (b') 

+ ^(aiia' (1) ) <g> (a' (2) * b) + ^(6i|b' (1) ) <8> (a* b' (2) ) 

(a') (b') 

+ 1 ® (at i (a' * b) + 6i 1 (a * b') + (a^i) 1 (a' * b')). 

On the other hand, by Eq. ([TBI . Eq. (I2ip and the induction hypothesis, we have 

A(a*b) = A(ai|(a'* b) + 6i|(o* b') + (ai6i)|(a'* b)) 

= 1 ® (ai| (a' * b)) + a x | A(a' * b) + 1 ® (6i 1 (o * b')) 

+6i 1 A(a * b') + 1 <g> ((aifei) 1 (a' * b)) + (aA) 1 A(a' * b) 
= 1® (ai|(a'*b)) +oi|(A(a')*A(b)) +1® (&i|(a*b')) 

+&! 1 (A(a) * A(b')) + 1 ® ((aifti) 1 (a' * b)) + (aA) 1 (A(a') * A(b')) 
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Then by Eq. fl20|) . the right hand side is 

fl l ■ (E ® a (2)) * (! ® b + 1 &'(!)) ® b'( 2 ))) 

(a') (b') 

+ 61 i ((1 ® a + 1 (a' (1) <g> a' (2) )) * ^(1) ® k(2)) 

(a') (b') 

+ (ai6i) 1 ( J]( a (i) ® a (2)) * 5Z( b (D ® b (2))) 

(a') (V) 

+ 1 ® (oii(a ; * b) + 61 1 (a* b') + (ai&i)i(a' * b')). 

This agrees with A(o) * A(b). 

Therefore %a is a bialgebra. By the definition of * and A, "Ha is connected filtered. 
Then !Ka is automatically a Hopf algebra by [21] Proposition 5.3], for example. □ 

We next give a differential version of quasi-shuffle Hopf algebras. 

Theorem 3.5. Let (A,d) be a commutative differential algebra. Extending d to 'Ka = 
®k>oAl k by defining, for a :— ai | ■ • ■ fa k G A^ k , 

k 

(22) d(a) = 5^a<,i I • ■ • |a i>fc , 

i=l 

a = I % ' J ^ 
lJ \ d(aj), J = i. 

T/ien (IHU, d) a differential Hopf algebra. 
Proof. We have 

k 

Aod(a) = ^ A(a ijl I • • • |a i>fc ) 

1=1 

= X^*' 1 1 l a i,j)®( a M+l I |Oi,fc)- 

i=l j=0 

Here we have use the convention that if j = 0, then a^i | • ■ • |a^j = 1 G k and if j = k, 
then a id+1 1 |a i;fe = 1 G k. 

With the same convention, we also have 



d o A(o) = d(y](ai I |a>)®(oj+i I |a fc ) 

fc 

= 2^ ( d ( a i I l%)®( a i+i I l«fc) + («i I • ■ ■ ia>j)®d(a j+ i | |a fc )) 



j=0 i=l 
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k 



i=j+l 

A o d(a). 



□ 



4. Renormalization of multiple zeta values and the differential 

structure 

We now apply the general setup in §[3] to the study of multiple zeta values. In order to 
apply the Differential Algebraic Birkhoff Decomposition in Theorem 13.21 we need to con- 
struct a commutative connected differential Hopf algebra (!K, d), a commutative differential 
Rota-Baxter algebra (R,P,d) of weight — f and an algebra homomorphism (p : "K —* R. 
We will provide these in § 14. f I which is modified with a differential twist from [29] for which 
we refer the reader for some of the details. We then use this decomposition to study the 
regularized and renormalized MZVs. 

4.1. Differential Algebraic Birkhoff Decomposition for multiple zeta values. 

4.1 .1. Directional regularized multiple zeta values. As mentioned in §[2], the expression 

I 



(23) C(S) = E 



Si St. 

n • • • flu 

n 1 >--->n k >0 1 K 

makes sense for integers s\, ■ ■ ■ , s& only when Sj > and si > 1. For other integers, we 
call the expression formal multiple zeta values since they are only formal expressions 
without a real meaning. In order to make sense of such formal expressions, we define the 
directional regularized multiple zeta values 

. prune _ n k r k e 

< 24 > *<GM== E „;,...'„;. 

ni>— >rife>0 1 K 

where r« are positive real numbers that are introduced so that the space spanned by the 
directional regularized multiple zeta values is closed under multiplication. 

Z([l];e) can be defined recursively as follows. Consider the following set of functions. 

For (s, r) G Z x M >0 and (e, x) G C x R with Re(e) < 0, define 



f^ s ^(e,x):= f([ s r ];e,x) = —. 
For vectors s = (sx, • • • , s&) G Z fc and f = (ri, • • • , r^) G (M>o) fc , define 

. , P (ni+x) rie . . . p (rik+x) r^e 

Then Z( [1] ; e, x) is also given by the recursive definition 

Z([;];e;*)=Q(/^(e;*)) 
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where Q is the summation operator [T7l HH] 

(26) Q(/)(x) = ^/(x + n), 

n>l 

and, for s = (si, ■ • ■ , s*.) G Z fc and r = (r 1; ■ • ■ ,r fc ) G (M >0 ) fc , 

Z( [J] ; e, x) = Q (f Sk (a, x)Z{ [£ [ [ [ ' ^ ];e,x)),k > 2. 

These are related to the multiple Lerch functions [TDl [T7| which are generalizations of the 
multiple polylogarithms 

2 ni . . . z nk 

\As x ,-,s k (z Xy - ■■ ,Z k ) = Yl 3k 

111 ' ' ' I ii, 

ni>"-nfc>0 K 

We then have 

Z([%e,0) = Z([^e). 



4.1.2. The (differential) Hopf algebra for regularized multiple zeta values. We next construct 
a Hopf algebra from the regularized MZVs to capture the algebra properties of these values. 
We consider the commutative semigroup 

(27) I (s,r)G ZxR >0 } 



with the multiplication 



■ r + r 

We similarly define the commutative semigroups 

(28) 9Jt+ = {/^ I (s,r) G Z >0 x M >0 }, OH = {/^ | (s,r) G Z <0 x R >0 }. 

For each of these semigroups 91, let = k 91 be the semigroup algebra. By Theorem 13 A\ 

^:=^(^)l" = ^k9T n , 

n>0 n>0 

with the quasi-shufne product and deconcatenation coproduct, is a connected filtered Hopf 
algebra. Further, for 9t = 9JT or 9Jt„, the operator 

diAn^Am, [ s r ] =r[ s - r 1 } 



is a differential operator. Thus by Theorem 13. 5[ 'K^ is a differential Hopf algebra. 
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4.1.3. The Laurent series for directional regularized multiple zeta values. Let Cj-fe}}^ -1 ] 
be the algebra of convergent Laurent series, regarded as a subalgebra of the algebra of 
(germs of) complex valued functions meromorphic in a neighborhood of £ = 0. Since 
ln(— e) is transcendental over C{{e}}[£] by [291 Lemma 3.1], we have 

the polynomial algebra with the variable T = — ln(— e) and with coefficients in Cjle}}^" 1 ] 
which embeds into the ring of Laurent series C[T][e _1 ,e]] with coefficients in C[T] by the 
remark after Lemma 3.2 in [29J. Thus we obtain an injective algebra homomorphism 

(29) uiC^Ke-^n-Cmie-Se]]. 
Then with the decomposition 

C[e-\e}\ = e- 1 C[e- 1 ]®C[[e]] 

of C^ 1 ,^]] into a direct sum of subalgebras, we have the vector space direct sum decom- 
position 

C[T][s-\e}} = C[T\[e-\e]] = e^pKe" 1 ] © C[T}[[e}} 

of the subalgebras e -1 C[:T] [e _1 ] and C[T][[e]]. Thus Cp1[e~\e]] is a Rota-Baxter algebra 
with the Rota-Baxter operator P to be the projection to e -1 C[7 1 ][e -1 ]: 



\n>N J k<-l 



This can also be directly verified as with the case of C[e 1 ,e]} (see [121 US]). 
The following facts are easy to verify. 

Proposition 4.1. (1) Define the operator 

d:C[T][e-\e\]^C[T][e-\e]\ 

by d(e h ) = ke k ~ x and d(T n ) = nT n ~ 1 /e. Then d is a differential operator on 
C[T][s-\e}]. 

(2) We have 

k 

(30) d(Z( [!] ; e, x)) = £ r t Z( ;s,x), 

i=i 

where e*j G Z fe is the i-th unit vector. 

(3) (C[e _1 , e}], P, d) is a commutative differential Rota-Baxter algebra in the sense of 
Definition \3.1\ 

Remark 4.2. We note that the differential operator d and the Rota-Baxter operator P do 
not commute in general. For example, d(P(T)) = <9(0) = ; but P(d(T)) = P(l/e) = 1/s. 



By [291 Theorem 3.3], 



*<[?);')= E 



gnirie ^n^r^e . . . e n k r k e 



ni>--->n k >0 



si So s k 
n l n 2 ' ' ' n k 
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has a Laurent series expansion in C[ln(-e)]{{e}}[e- 1 ] for s G Z fe , f G Z^ and a Laurent 
series expansion in C{{e}}[£~ 1 ] for s G Z< , f G Z> . Combining this with Eq. (129j) . we 
obtained an algebra homomorphism 

Z : - C[T] [e- 1 , e]], ^ ) = u{Z{ [1] ; e)) 

which restricts to a differential algebra homomorphism 

Z : % m - -> C[e _1 ,e]]. 
Then by the Differential Algebraic Birkhoff Decomposition in Theorem 13.21 we have 
Corollary 4.3. There is the unique decomposition 

Z = ZZ l * Z + 

where the map Z + : "Km ~~ * C[T][[e]] zs an algebra homomorphism which restricts to a 
differential algebra homomorphism Z + : "H^n- ~ * C[[e]]. 

By Theorem 3.8 in [29J, 

(31) Z + ([-«]; £ ) = (id-P)(z([« ; «]; £ )-P(Z([3; £ ))Z([«]; £ )). 

Definition 4.4. For s = (si, ■ ■ • , s*.) G Z fc and r = (r 1; ■ ■ • , r&) G R> > define the renor- 
malized directional multiple zeta value (MZV) 

(32) CO=limZ + (&4 

The definition makes sense because of Corollary 14.31 Furthermore, since Z + is an algebra 
homomorphism, we have the quasi-shuffle relation 

c([3x([?])=<([a*[?]) 

meaning if 

is the quasi-shuffle product in 'Kyi defined in Eq. ( TTBl - (ITT)) , then 

c([JlK([j!])= Ec([?D- 




Definition 4.5. For s G Z> U Z< , define 

(33) C>1=£mC([ |s1 s ;j), 

where, for s = (si, • • • , Sfc) and 5 G R>o, we denote \s\ = • • • , |sfc|) |s| + 5 = 
(|si| + 5, • • • , |sjfc| + 8). These are called the renormalized MZVs of the multiple zeta 
function ((ui, • • • , Uk) at s. 
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We show in [29l Theorem4.2] that our renormalized MZVs are well-defined and are com- 
patible with the known MZVs defined by either convergence, analytic continuation [46] or 
the Ihara-Kaneko-Zagier regularization [33]. It is also proved that the renormalized MZVs 
satisfy the quasi-shuffle relation. 

We next show that our examples in Section [2] follows from the above general set up. Let 
Sfc denote the symmetric group on k letters. For a G E& and r — (ri, ■ • • , r^), denote 
cr{r) = (r CT( i), ■ ■ • ,r a[k) ) and f(r)^ = £ ffeSfc f{<r{r))- 

Proposition 4.6. ([29], Proposition 4.9]) For any k>l, C([°^])^ S ^ i s independent of the 
choice of f e R> and 

Remark 4.7. Taking k = 2 and r = (1, 1) in the above proposition, we have ((0, 0) = 
C([i'i])- Then by Eq. (I3T1) and Eq. (1321) . we obtain the renormalized value for C(0, 0) 
discussed in Section |2j 

4.2. Differential structures on renormalized multiple zeta values. We now consider 
some differential properties of the renormalized MZVs. By Eq. ([SI), 

*<C]i«> = ~+E«-o^. 

i=0 



So 

(34) Z + ([° r ];e) = J2a 

8=0 

is the generating function for C( s )> s £ ^<o- We next generalize this to multiple variables. 
Theorem 4.8. 

(35) = E E <([-<"••,• • fe >])^-^, 

n>0 iiH \-ik=n 

where the sum runs over ordered partitions of n. 
Proof. Let 

n>0 

We find a formula for a n . Of course a n is the constant term of 

d n ~ 
de r 



But since Z + is a differential algebra homomorphism by Corollary 14.31 we have 

-£m;;h=e-^+([*V'h 



de 

i=l 
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where is the i-ih unit vector of length k. By an inductive argument, we have in general 

z { :\^e)= £ ( jii : iii M-ft([ ? -^f'' ,t, ]i e ) 

iiH h«fc=n 

where the sum runs over ordered partitions of n. Thus we have 

r(n) / rs" 



n nj (C-];o) 

iiH hifc=n 



£ ( h ".K-^z + ([- l V' ,, "]i'» 

-+ife=n 



ii,— * 

HH h'tfc=n 

Therefore we have the generating function 

n>0 iiH h«fe=™ 

E E cd'-'V''*']) 



M3) = E E ( h ; A K"-*c(['- ( V' w ]) 



2X ifc 

S- (ii, • • • ,i fc )l r j r fc c n 



Zi! 2z~! 
n>0iiH M fe =re 

=E E c(['-<V-«])^-^ 

n>0 iiH Mfc=n 

This proves Eq. ( 1551) . □ 

We also consider the regularized series before the renormalization. In the one variable 
case, we have 



1 oo 

For the two variable case, we have 



re ' — ' z! 

i=0 



2(0^) = P »+E"4 



n>0 



where Po is the negative power part. Let us find a formula for a n . First by n-th derivative, 
we have a n =constant term of But 



i=0 



and, by [29l Proposition 4.8], the constant term of Z([ ^^"jje) is 
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So 



E(D^rc([l:; n ]) 

i=0 
n 

EO-i'-r'cdt!;"]) 



n 

v i z + 1 ri 



„n+l 



i=0 
,n+l 



C(-n-l) 



i=Q 



n 
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